Towards higher precision and operational use of optical homodyne tomograms 
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We present the results of an operational use of experimentally measured optical tomograms to 
determine state characteristics (purity) avoiding any reconstruction of quasiprobabilities. We also 
develop a natural way how to estimate the errors (including both statistical and systematic ones) 
by an analysis of the experimental data themselves. Precision of the experiment can be increased 
by postselecting the data with minimal (systematic) errors. We demonstrate those techniques by 
considering coherent and photon-added coherent states measured via the time-domain improved 
homodyne detection. The operational use and precision of the data allowed us to check for the 
first time purity-dependent uncertainty relations and uncertainty relations for Shannon and Renyi 
entropies. 
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INTRODUCTION 



A measurement plays a vital role in the study of quan- 
tum physics. Optical homodyne tomography is merely 
one among a variety of measurement techniques, how- 
ever, its importance and effectiveness can scarcely be 
overestimated. The conventional optical homodyne to- 
mography of one-mode continuous-variable states takes 
its origin from the papers [l|-[3] and is instructively de- 
scribed in a series of books and reviews (see, e.g., 0- 
6]). The original goal of the optical homodyne tomog- 
raphy was to infer the quantum state of light identi- 
fied with the density operator p or the Wigner function 
W(q,p) p}, say. In fact, any faithfully reconstructed 
quasi-probability contains a complete information about 
the state and can then be used to calculate any charac- 
teristics of the state, for example, its purity. Unfortu- 
nately, no reconstruction procedure is perfect and, what 
is more unpleasant, the original errors of experimental 
data can grow during the reconstruction. It is generally 
accepted that the higher precision of the measurement, 
the more comprehensive information is provided and the 
more sophisticated phenomena can be observed. The pre- 
cision is thought to be increased merely by increasing the 
number of experimental runs (enlarging an ensemble of 
identically prepared states). Clearly, such an approach 
leads to a reduction of statistical errors but can hardly 
cope with systematic ones (related with the experiment 
itself). On the other hand, quantum tomography is a 
quantitative technique only if we can evaluate the overall 
errors presented in the experimental data. The previous 
approaches do not give a direct solution of this problem: 
the pattern-function reconstruction can provide the sta- 
tistical errors only, whereas the maximal likelihood ap- 
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proach to evaluation of the errors resorts to a bootstrap 
method whose result cannot be totally relied on and is 
time-consuming [6]. In this paper, we propose and apply 
in practice a straightforward method to estimate both 
statistical and systematic errors, thus, making a step to- 
ward a higher precision of homodyne experiments. 

Our approach is based on the initial experimental data 
and corresponding histograms of quadrature values. The 
histograms are nothing else but an estimation of the 
quantum tomogram w(X,6). Being a measurable char- 
acteristic of the state and describing the quantum state 
thoroughly, the tomogram is not only a powerful tool to 
reconstruct quasi-distributions but can be used solely on 
an equal footing (see the tomographic-probability rep- 
resentation of quantum mechanics [8|, |9(). Moreover, as 
we show in this paper, the tomographic approach can 
be used to estimate the errors of histograms and, what is 
more important, to calculate directly state characteristics 
(e.g., the purity) and their errors. 

Although different kinds of photon states can be an- 
alyzed by the optical homodyne tomography, we focus 
our attention on photon-added states [10l4l2l |. whose 
experimental detection |13l - fl7| and nonclassical behav- 
ior [13, [l9j were demonstrated recently. Moreover, the 
advanced techniques of photon addition and photon sub- 
traction enabled us to perform a direct probe of the com- 
mutation relation between photon creation and annihila- 
tion operators (23423 as well as accomplish a noiseless 
amplification [23| . 

We use coherent (classical-like) and single-photon 
added coherent (non-classical) states to achieve another 
goal of our paper, namely, to analyze the accuracy with 
which the known so far fundamental quantum relations 
are fulfilled. Such relations include, for example, the 
Heisenberg inequality 2A\ and its purity-dependent ver- 
sion 12511 as well as the state-extended uncertainty rela- 
tions [2a, [23| and the uncertainty relations for Shannon 
and Renyi entropies |28l - l30| . It was shown theoretically 
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FIG. 1: (Color online) Schematic of the experiment for the 
generation of SPACSs. An UV pulse pumps a nonlinear crys- 
tal to produce stimulated parametric down-conversion (PDC) 
in the mode of a seed coherent state. Detection of a single pho- 
ton in the conjugated idler mode by a single-photon counting 
module (SPCM) heralds the successful generation of a SPACS 
in the signal mode, and triggers its homodyne detection. This 
is performed by mixing the signal state with a coherent local 
oscillator (LO) pulse on a 50% beam-splitter and measuring 
the difference photocurrent produced from two photodiodes 
at its outputs. 



in the papers [31r436l | how to check all these inequalities 
by means of the optical homodyne tomography. In this 
paper, we present the first experimental results for some 
of them. Needless to say that the accuracy of tomo- 
graphic data plays the major role in this case. However, 
a fulfillment of the Heisenberg uncertainty relation does 
not mean that the quantum mechanics in its conventional 
form is valid (see, e.g., [33, [38[) and opens a possibility of 
goin g be yond the conventional quantum mechanics (see, 
e.g., |39|,|4(|). The violation of the conventional quantum 
mechanics (if any) could be detected by the violation of 
quantum inequalities for highest moments. In principle, 
all the highest moments can measured via homodyne de- 
tector as well [41| and the experimental check of those 
inequalities is to be discussed elsewhere. 

The paper is organized as follows. 

In Sec.|lTl the optical homodyne tomography is shortly 
reviewed with the emphasis on coherent and photon- 
added coherent states. Also, an optimal estimation of 
the quantum tomogram is developed and the influence of 
detection imperfection is discussed. Sec. [Hi] is devoted to 
the estimation of errors and a brief analysis of the reasons 
of systematic errors. In Sec. IIV1 we describe how to deal 
with the experimental data in an operational way by cal- 
culating the purity of the state and checking uncertainty 
relations mentioned above. Conclusions and prospectives 
are given in Sec. [V] 



II. OPTICAL HOMODYNE TOMOGRAPHY 

The basic idea of the homodyne tomography is to mea- 
sure the quadrature operator Xg — QcosO + PsinO, 
where Q and P play the role of position and momentum 



such that [Q, P] = ih and 9 £ [0, 2tt] is a phase of a strong 
coherent light also called the local oscillator (LO). Note 
that Q and P have the same units and h is a constant, 
specified during the calibration procedure. Fixing the 
LO phase 9, one can get access to the probability density 
distribution (tomogram) w(X,9) — (Xg\p\Xg), where 
Xg\Xg) — X\Xg). If the tomographic values w(X, 9) are 
specified for all the points X £ (— oo, -t-oo) and 9 £ [0,7r), 
then such an ideal tomogram contains the complete in- 
formation about a quantum state. 

In this section, we show how to estimate the tomogram 
of a coherent state and a single photon added coherent 
state (SPACS) in the experiment. Before we move on to 
the description of the experiment we briefly discuss the 
states under investigation. 



A. Coherent and SPAC states 

Coherent state |cv) is an eigenstate of the photon an- 
nihilation operator a = (Q + iP)/\/2h, viz., a\a) = a\a), 
where a £ C. The coherent state \a) is determined by 
the following tomogram: 



W| a >(X,0) = ^i=exp{- 



■y? — v2 (Rea cos 9+Ima sin ( 



(1) 

The SPACS is defined as a t |a)/ v /l+ |a| 2 and its to- 
mographic representation reads (see, e.g., |42j |) 

wu\ a) (X,6) = [v^fi(l + \af 

-y= — i= (Rea cos 9 + Ima sin t 
+(Rea sin 9 — Ima cos 9) 2 > 
xexpj- 77= - V2(Rea cos 9 + Ima sin 9) j. (2) 

It is not hard to see that in the limit |a| — > oo formula 
© reduces to ([1]), i.e. the SPACS behaves as a coherent 
state. A transition from a purely quantum behavior of 
the SPACS (a = 0) to a classical-like one (\a\ ^> 1) was 
also observed experimentally [13j. 



B. Experimental setup 

SPACSs are generated by injecting a coherent state 
| a) into the signal mode of an optical parametric am- 
plifier and exploiting the stimulated emission of a single 
down-converted photon into the same mode. Successful 
SPACS generation takes place upon detection of a single 
photon in the idler mode of the amplifier. Quadrature 
data are then acquired by a time-domain balanced ho- 
modyne detector [15l 1 161 ] triggered by such idler counts. 
A schematic of the setup, described in detail in [13l.ll4J. 
is presented in Fig. [1] 



(a) 



(b) 



(ci 



1.0 

as 
0.6 
0.4 

0.2 



ill.: 



ulfellUll 




1.0 






0.8 


— 




0.6 


i- 


1 


0.4 


r 




0.2 


— r-r-i-rrTI. 


. ,h-u_. 



1 2 



FIG. 2: (Color online) Histogram h(X, re) of a SPACS for different bin widths: (a) b = 0.025, great statistical errors; (b) 
b — 0.075; (c) 6 = 0.15, larger bin widths can cause undersampling. 



Acquisition of the quadrature data from the homodyne 
detector is accomplished by means of a digital oscillo- 
scope, producing a sequence of N = 5321 quadrature val- 
ues X for each fixed LO phase. Calibration of X- values 
is accomplished by measuring vacuum fluctuations when 
the signal is blocked. In this case, (X) = and the vari- 
ance axx = ((X — (X)) 2 ) = ft/2. Thus, a choice of h is 
rather arbitrary and we use H = 4. Once X is calibrated, 
a state under investigation is characterized by a collec- 
tion of points {Xi, 9j}, where i = 1,...,N. The phases 
9j are adjusted by the piezoelectric transducer. 



C. Tomogram estimation 

The binned histogram h(X,9j) is known to be con- 
structed ambiguously because of many possibilities to 
choose the bin width b. If b — > 0, then the histogram 
is merely a sum of delta functions S(X — Xi). In fact, for 
relatively small bin widths no statistical confidence can 
be achieved. Conversely, if 6 — > oo, then the histogram 
transforms into a flat distribution over the range of X , 
with this uniform distribution tending to zero. In this 
case, no useful physical information can be extracted. 
Needless to say that none of these two extremal types 
of the histogram reflects the behavior of the function 
w(X, 9) predicted by the theory. 

Let us now derive an optimal bin width b for pur- 
poses of the optical homodyne tomography. To begin 
with, the histogram value h(Xi,9j) at the point Xi = bi, 
i £ Z, equals Ni/Nb, where Ni is the number of measured 
quadrature values falling into the i-th bin [Xi,X{ + b) 
and N is the number of all quadrature values. The 
statistical error of h(Xi,9j) originates from Ni whose 
error is ^/Wl because the measurement process is as- 
sumed to be Poissonian. For a fixed N we naturally 
have Ni oc b (if b is not too large), then the statisti- 
cal error of h(Xi,9j) is 5h stajt — y/h(Xi,6j)/Nb, For 
relatively large bin widths the statistical error is negli- 
gible, however, the effect of undersampling the quadra- 
ture distribution comes into play [43j. The main idea 
is that the theoretical tomogram w(X,8j) exhibits os- 



cillating behavior with respect to X, with the scale of 
oscillations being ~ 7r/v2d, where d is the number of 
Fock states significantly contributing to the state un- 
der investigation. The experimental histogram h(X,9j) 
should reflect those oscillations rather than conceal them. 
Then, the error of undersampling for the histogram value 
h(Xi,9j) can be evaluated as 5h nn d — h(Xi,6j)bv2d/n. 
The resulting error <5/i s tat + Sh un( { takes minimal value 
if b = b opt = [n/iV^hiX^d^Nd] 1 / 3 . Note that b has 
the same functional dependence ex 1/ y/N as the Scott's 
choice b = 3.5<r / \7~N , where a is the standard deviation 
of X [44J. Note also that the optimal bin width should 
increase for lower values h(Xi, 9j), for instance at the end 
of the distribution tails. For practical purposes the al- 
ternating bin widths are, however, not very convenient 
since they complicate data processing. 



We plot some examples of histograms for different val- 
ues of the bin width in Fig. [5] For our further anal- 
ysis, we choose b — 0.075 which is close to the av- 
erage optimal value b opt ~ 0.06 (we put h(Xi,9j) f=s 
l/v2vr, N = 5321, d ~ 1, and scale 6 op t by a factor 
\/h = -4s). In our case, this bin width is also close to 

(max X — mm X) / i/N « 0.055 known as the square-root 
choice. For the normal distribution (coherent state) the 
Scott's choice gives b = 0.14 and the Sturges' formula 
results in [log 2 N + 1~\ w 13 bins (b — 0.3). We choose 
b = 0.075 to guarantee the statistical confidence and pre- 
vent the data from undersampling. The latter fact is 
important to observe the cases when the theoretical func- 
tion w(X,9) tends to zero in the middle of the range of 
X as it takes place, e.g., for an ideal SPACS (see formula 
([2]) with real a and 9 = or if). 



To illustrate the estimated tomograms of a coherent 
state and a SPACS, we present a series of historgrams 
h(X,9j) constructed on the basis of the experimental 
data with eleven LO phases {0j]4jLi within the region 
[0,tt] (see Fig. ED. 




FIG. 3: (Color online) Typical histograms h(X,8j) of a coherent state (left) and a SPACS (right). 



D. Detection imperfection 

Let us be reminded that for a real a and the LO phase 
9 = the tomogram fl2J of a SPACS is tu at | a > (X, 0) oc 
(X/\fh—a/\/2) 2 exp[— (Xj\fh~ y/2a) 2 ] and takes on zero 
value if X = a^/h/2 (X = —a^h/2 it 9 = n). However, 
one can hardly observe such property in Figs. [2] and 
This is caused by the fact the detection efficiency rj < 1. 
The detection efficiency comprises all kinds of losses in- 
cluding the finite efficiency of photodetectors. Due to 
the imperfect detection, the measured histograms are 
smoothed and there is no zero point X = a^h/2 any- 
more. In fact, one actually measures not the prepared 
state but its convolution with a vacuum (that impinges a 
fictitious beamsplitter with transmittivity r\ in front of an 
ideal quadrature detector) . In terms of the Wigner fac- 
tion, the measurable state W det (q,p) is connected with 
the originally prepared state W(q',p') by the following 
relation: 



W aet (q,p) 



1 



7T(l-7?) 



dq' dp' W(q',p') 



x exp 



(<z-vW) 2 + Cp-vV) 2 



1-7J 



.(3) 



It is not hard to see that a coherent state |a) transforms 
into the coherent state \y/rja) under convolution ©. On 
the other hand, a SPACS remains no longer a SPACS and 
the measurable state is given by the following Wigner 



function: 

<t°| t Q) (?,p) = Ni + |a| 2 )]- 1 



x{i+ 2 #-^H +(?-sH - 1 ]} 



x exp 



2?7Rea: 



P- 



27/Ima 



(4) 



Then, for a SPACS with real a, the theoretical prediction 
of the measurable quadrature distribution is 

-, -l 



..dot 



#f a) (X,0)= Vnh(l + a') 



i-^ 2 ^- 2 ^)' 



exp 



TTTr-VV* 



•JTi 



(5) 

which has no zeros and correctly describes the experi- 
mental histograms in Figs. [5] and [3J 

It is worth noting that the purity of the state can reveal 
the detection imperfection. Although a coherent state 
remains pure in transformation ([3]), a SPACS does not. 
Indeed, the purity of the detectable SPACS reads 



det 

Matia) 



= 2tt 



dqdp W£?Uq,p) 



= 1 - 



21,(1-77) 



(i + H 2 ) 2 

(6) 
and is less than 1 whenever 0<?7<l(if?7 = 0, then the 
vacuum noise is only detected). 

In what follows, we will concentrate on the accuracy 
of the experimental histograms h dct (X, 9) and theoreti- 
cal tomograms w dct ( X, 9) . Thus, we will operate with 



the "detectable" state (not the originally prepared one). 
Further, we will omit the superscript dct wherever it is 
clear from the context. In fact, deconvolution of formula 
([5]) is known to be difficult to perform with experimen- 
tally given quasiprobabilities [J] and this is beyond the 
scope of present paper. We can refer the interested reader 
to the paper [45| , where a similar deconvolution problem 
is solved, namely, an extraction of the originally prepared 
microwave quantum state from a noisy output of a linear 
amplifier is considered. 



III. ACCURACY OF OPTICAL HOMODYNE 
TOMOGRAMS 



TABLE I: Shifts x = (X n ) + (X ) of the histogram h(-X,0) 
with respect to h(X, n) for detected coherent and SPAC states 
of different intensities. The amplitude ^/rja of the detected 
coherent state is evaluated by the experimentally measured 
value ((AT horcnt > - <X§ ohcrcnt ))/2. 





Data 


Detected amplitude ^/rj 


a 


set 


0.64 


0.82 


1.25 


1.73 


#1 


coherent 


0.14 


0.15 


-0.15 


-0.17 


SPACS 


0.16 


0.20 


-0.08 


-0.10 


#2 


coherent 


0.26 


-0.21 


0.02 


0.23 


SPACS 


0.26 


-0.14 


0.05 


0.27 


#3 


coherent 


0.03 


-0.12 


0.003 


0.36 


SPACS 


0.09 


-0.07 


0.07 


0.40 



Further progress of applied quantum information tech- 
nologies and fundamental experiments depends greatly 
on the accuracy of measurement data. In optical ho- 
modyne detection of radiation held, one usually restricts 
oneself by the initial calibration of the detector outcomes. 
Namely, blocking photons of the signal mode results in 
the vacuum state, whose quadrature distribution is to 
be centered at point X — and have the dispersion 
(X 2 ) = tl/2 for any phase of the local oscillator. How- 
ever, in practice, a drift of the scheme parameters or an 
extra noise can occur during the experiment. In view 
of this, for practical purposes it is extremely important 
to trace the adequacy of the data being collected either 
in real time or during postprocessing. Also, the method 
would be beneficial if it were based on the data them- 
selves without much additional information. In this sec- 
tion, we present and apply such a method. 

A true tomogram w(X, 9) is known to satisfy the rela- 
tion w(X,9) = w(—X } 9 + 7r). This fact was previously 
used to claim that the quadrature distribution for LO 
phases 6 £ [0,7r) determine a quantum state thoroughly. 
As a result, the phases out of this range were disregarded 
in experiments, although they naturally provide an effi- 
cient way to check the accuracy of the data. In what 
follows we show that one can efficiently use the peculiar 
property w(X, 9) — w(—X, 9 + n) to check whether the 
data are adequate (46|. Moreover, one can evaluate the 
accuracy of the histograms. 

For example, an imbalance of the optical scheme or 
photodetectors' efficiencies would result in values X 
shifted by some cc; m b. In this case, the distributions 
w(X,9) and w(—X,9 + %) as functions of variable X 
would be shifted with respect to each other by the magni- 
tude 2x; m b. In case of different photodetector efficiencies, 
r\\ and 77 2 , the shift x\ m \, oc (771 — 772)/, where / is the LO 
intensity. Analogous mismatch between tomograms can 
take place due to a low frequency electronic noise at the 
input of the digital scope, the shift alternating in time. 

Another reason of possible deviation of w(X, 6{) from 
w{— X, 62), where 9 2 is supposed to be equal to 9\ + n, 
can occur due to inaccuracy in the LO phase control. 
Especially clearly this type of data mismatch is seen for 
a coherent state \a), for which the distribution w{X, 9\) is 



shifted with respect to the distribution w(X, 62) by xse = 
V / 2^[Rea(cos 9\ + cos 62) + Ima(sin 9\ + sin 9 2 )) along X- 
axis. For #i, 2 = 9 ± 89/2, the shift xse is approximately 
equal to \/2h(J{ca sin 9 — Ima cos 9)89 + h(Reacos9 + 
Imasm9)89 2 /V2. 

In order to demonstrate the method above, we consider 
a mismatch between histograms h(X,ir) and h(—X,0), 
which should be coincident according to the theory. Typ- 
ical histograms of a SPACS are depicted in Fig. @] for 
three data sets corresponding to ^/rja = 0.64. 

One can readily notice the deviation of histograms 
h(X,%) and h(-X,0) for data sets #1 and #2. The 
shift between these histograms is evaluated as the dif- 
ference between mean values of the distributions (to be 
precise, the shift x = (X n ) — (—Xq) = {X„) + (X )). 
The experimentally determined shifts of the histogram 
h(—X, 0) with respect to h(X, ir) are summarized for co- 
herent states and SPACS states of different intensities in 
Table H 

While I a I is getting larger, one expects the error of fix- 
ing the LO phase 89 to get smaller (since the phase con- 
trol is based on observing an interference picture which 
becomes clearer for larger |a|). On the other hand, in 
case of real a and LO phase 9 = 0, the shift xse equals 
ha89 2 J 'v2 and can be non-monotonic with respect to a 
because of an additional factor. 

Let us now analyze how a mismatch between distribu- 
tions w(X,9) and w(—X,9 + 7r) affects the accuracy of 
the data and allows evaluating the experimental errors of 
some state characteristics. 

A natural characteristic, which shows the closeness 
of two probability distributions Pi(X) and p 2 {X), is 
the Bhattacharyya coefficient [47( defined as B = 
J y/pi(X)p 2 (X)dX . The Bhattacharyya coefficient B 
equals 1 if and only if distributions Pi(X) and p 2 (X) 
are identical. 

Let p\ be a state reconstructed from the homodync 
tomograms w(X,9), 9 <E [0, 7r), and p 2 be a state re- 
constructed from the tomograms w(X,9), 9 E [tt,2tt). 
Provided ideal tomograms the states p\ and P2 are iden- 
tical. Experimental data result in two different states, 




FIG. 4: (Color online) Histograms h(X, n) (blue solid lines) and h(—X, 0) (red dashed lines) of SPACS from first column of 
Table |T] (a) data set #1; (b) data set #2; (c) data set #3. 



the fidelity F(pi,p2) — Tr -J y/piP2 \fp\ between which 
indicates the accuracy of measured data and can be 
used as an estimate of the fidelity between the evalu- 
ated (reconstructed) state yO es t and the actual state p, 
i.e. F(p\ 1 p2) ~ F(p,p cst ). Important for us is the fact 
that F(pi,p2) satisfies the following relation 48]: 



F < min 

9G[0,7r] 



y/w(X,9)w(-X,8 + 7r)dX, (7) 



that is the fidelity is limited by the minimal Bhat- 
tacharyya coefficient Bg for the distributions Pi(X) = 
w{X,9) and p 2 (X) = w(-X, 6 + n). 

In principle, formula implies minimization of Bg 
over all experimentally accessible LO phases 9. In this 
research, we restrict ourselves by an illustration of the 
method of fidelity evaluation and present some values 
Bg = o calculated for the data from the first column of Ta- 
ble ED 98.70% and 98.32%, 96.20% and 95.59%, 99.67% 
and 99.26% for the coherent and SPAC states from data 
sets #1, #2, and #3, respectively. Here, we have cal- 
culated the integral (0 by replacing w(X, 9) — > h(X, 9) 
and using the trapezoid method [49j , with the error of cal- 
culation being --^b 3 -^ v / w{X,0)w(~X,Tr) < 0.004%. 
Once fidelity is evaluated, one can use this knowledge to 
evaluate the accuracy of other state characteristics (see, 
e.g., SHI). 

Given tomograms w(X, 9) for two regions of the LO 
phases 9 e [0,ir) and 9 £ [tt, 2tt), it is possible to eval- 
uate the error of the mean value of any physical quan- 
tity A. Indeed, AA = |Tr[(pi - pi)A\\, where p\ and 
P2 are defined as above. However, for some quanti- 
ties one does not have to reconstruct the states and 
can use tomograms directly. For instance, the moment 
{Xq) = J X n w(X,9)dX is determined with the experi- 
mental error A(A^) = / X n \w(X, 6)-w(-X,6 + n)\dX. 
For example, for the data set #3 from the first col- 
umn of Table U the second moment (q 2 ) = (Xg = o) 
equals (q 2 ) = 0.63 ± 0.04 for the coherent state and 
(q 2 ) = 1.23 ± 0.16 for the SPACS. For the first moments 
(q) the errors are merely the shifts 0.03 and 0.09, respec- 
tively. The error bars of those quantities are of the same 
order for other LO phases. 



To conclude this section, a relatively simple analysis 
of the homodyne tomographic data enables one to check 
their adequacy and evaluate their accuracy. As a re- 
sult, one can postselect and use further only those data 
that meet the desired accuracy. Moreover, a mismatch 
between tomographic data can indicate a reason and na- 
ture of extra noise. The latter fact opens up new vistas 
of the optical homodyne tomography in metrology. 



IV. OPERATIONAL USE OF THE 
TOMOGRAPHIC DATA 



In this section, we are going to reveal some relevant 
information about a quantum state just using the tomo- 
graphic data and circumventing a reconstruction of the 
density operator or the Wigner function. Also, we are 
checking if the data satisfy some theoretically predicted 
inequalities. In this section, explicit numerical values of 
quantities of interest are calculated for data set ^3 from 
the first column of Table U which exhibits relatively small 
systematic errors. 



A. Purity 

Purity p = Trp 2 is an important state characteristic 
which can set some limitations on the use of the state in 
applications. A conventional approach to determine the 
state purity from the optical tomogram is to reconstruct 
the density matrix or the Wigner function via some im- 
proved modifications of the inverse Radon transform 52 1 
or the maximum likelihood method [53J and then sub- 
stitute them in some integral relations to calculate the 
purity. Recently, the state purity has been also evalu- 
ated from quadratures' uncertainties J5J|. This method 
is easy to use but the evaluation gives a correct value 
only for Gaussian states. Here, we use the tomographic 
data directly and calculate the true purity without any 
intermediate reconstruction of the density operator or a 
quasiprobability distribution. Moreover, no assumption 
about the state being Gaussian is needed. 



The purity is known to be expressed through the opti- 
cal tomogram as follows |55| : 



1 



Mtho 



+00 



+00 



drr I I dXdYe^ x+Y ^ r 
2tt J J J -00 

x / d9w(X,6)w(-Y,9), 



(8) 



where the sequence of taking integrals is chosen for the 
easiest data processing. If the tomograms satisfied the 
relation w(X, 9) — w(—X, 9 + n), the calculated value of 
\x would be real. In fact, one would have 




d9w(X,9)w{-Y,9) 
d9 [w(X, 9)w(-Y, 9) + w(-X, 9)w{Y, 9)] (9) 



and, consequently, 



1 p+00 p p+00 

H = - drr dX dY cos[(X + Y)r] 



x / d9w(X,9)w{-Y,9). 

Jo 



(10) 



The obtained formula is beneficial when the homodync 
data are acquired only for the LO phases in the range 
[0,7r] (although it is impossible to evaluate the accuracy 
of /i then). As we already know, in practice the require- 
ment w(X, 9) = w(— X, 9 + 7r) is not precisely met. Then 
the imaginary part of expression ([SJ can serve as the error 
bar of the purity. It can be also calculated as follows: 

A M = (Trp 2 - TrpD/2 

1 p+00 p p+00 

= — rdr dXdY cos[{X + Y)r] 

2tt J JJ_ 00 

i9[w(X, 9)w(-Y, 9) - w(X, 9 + n)w(-Y, 9 + n)}. 

(11) 

Given experimental histograms h(X,9j), we first 
calculate the sum ^Ylj=i [h{Xi,9j)h{—Y k ,9j) + 
h(Xi,9j-\-i)h(—Yk,9j+i)](9j+i — 9j) for any pair of bin 
coordinates (Xj,Yfe), i.e. the evaluation of the function 
P(X,Y) = f* d6w(X,9)w(-Y,9) via the trapezoid 
method. The error of this evaluation is roughly 
equal to -±(9 j+1 - 9 j ) 3 JL[w(X i ,9)w(-Y k ,9)} < 



for the states 



Xf 



Y k 2 ) 



< 



0.003 exp(-X| - Y 2 ) 



in question. Calculation of the in- 
tegral J(r) = JJ^dXdY cos[(X + Y)r]P(X,Y) 
is substituted by the calculation of the sum 
J2x t ,Y k b2cos i(X l + Y k )r]P{Xi,Y k ) for any fixed r. 
This evaluation contains two types of errors: the first 
one originates from the error of the function P(Xi,Y k ) 

and equals 0.01e~ r / 2 , and the second one is due to evalu- 
ation of the integral JJ_ dXdY by the sum J2x y an( ^ 



FIG. 5: (Color online) Calculation of purity of the SPACS. 
The purity is the area under the curve divided by 2ir. Func- 
tion J(r) deviates from for r > 8 due to calculational errors. 



< 



equals -^6 4 (^ + ^) cos[(X l + Y k )r]P(X t ,Y k ) 

10 _5 (2 + r 2 ). Evaluation of the function rJ(r) for the 
SPACS is presented in Fig. [5j Deviation of J(r) from 
for values r > 8 is to be assigned to the second type of the 
error. Finally, the purity parameter (|10[) is calculated via 
integrating the function rJ(r) in the range [0,-R], where 
the upper limit R is chosen in such way that the integral 
/i(i?) := L rJ(r)dr is saturated and does not depend 
of R. The error of calculating p(R) can be evaluated as 
A ca i c /i(i?) < 0.01 f+°°re- r2 / 2 dr+10- 5 f* ' r(2 + r 2 )dr a 
0.01 + 3 • 10- 6 i? 4 . We choose R = 8 and obtain /i = 1.00 
for the coherent state and fj, = 0.83 for the SPACS, the 
error of calculation being A ca i c /i < 0.02. Similarly, we 
use formula (jlip to calculate the error A/z originating 
from the inaccuracy of the original data. The direct 
calculation yields 0.035 for the coherent state and 0.039 
for the SPACS, which is slightly greater than the error of 
calculation A ca j c //. To resume, we obtain fi = 1.00±0.04 
and n = 0.83 ± 0.04 for the coherent state and the 
SPACS, respectively (the overall error is estimated as 
[(A calc ^) 2 + (Afi) 2 } 1 / 2 ). It is worth noting that the 
obtained purity fj, = 0.83 of the detected SPACS exactly 
coincides with what is predicted by Eq. © if nominal 
values 77 = 0.6 and a — 0.83 are used. 

We note that the errors are easily and naturally es- 
timated in our approach in contrast to the approaches 
based on the density operator reconstruction that in- 
volves a rather time-consuming bootstrap method for 
evaluation of the errors by the maximal likelihood tech- 
nique. Note that a calculation of the purity for a given 
density operator also results in additional calculational 



B. Fidelity 

Usually, an experiment is aimed at producing a spe- 
cific pure quantum state, \ip) say. Experimentally de- 
termined quadrature distributions w e xp(A, 9) allow cal- 
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culating the fidelity F 2 — (ijj\p\ip), where p is an ac- 
tually detected state. Similarly to formula (fTDj) . we 
have [5l: &W) = ± J Q + °° r dr ff+™ dXdY co S [(X + 
Y)r] Jq d0w^(X,6)w exp (—Y,6), where the analytical 
function w^{X,0) — |(X6/|-0)| 2 is easily computed 
through the desired state \ip). For instance, if 
\ip) is a superposition of a finite number of Fock 
states \n), then an explicit formula for w^{X,9) 
is found, e.g., in Ref. [4y|. Since the func- 

tion Wij, (X, 9) is known precisely, the error of the 
quantity F 2 equals ^ L r drff_ 00 dXdY cos[(X + 
Y)r] /„" dOw^X, 6)[wl p (-Y, 9) - w c Z(Y, 9 + n)}. 



C. Experimental check of uncertainty relations 

1. Heisenberg inequality 

Since the main difference between two histograms 
h(X,8) and h(—X,9 + n) is essentially the shift, the 
variances <Jx g x e — (A|) — (Xg) 2 differ not so severely 
as the second moments. In fact, in our case we have 
Aa qq — 0.004 for the coherent state and 0.013 for the 
SPACS. 

In this subsection, we are going to check if the Heisen- 
berg uncertainty relation a qq a pp > h 2 /A holds true and 
what is the extent to which it is fulfilled. Certainly, 
thanks to the initial calibration of the apparatus by the 
vacuum state, we can adjust H = h and check the in- 
equality for any other states but not the vacuum itself. 
The errors of determining second moments are evaluated 
as described in Sec. IIIII 

For the coherent state we have a qq a pp = 0.0612 ± 
0.0014, which coincides with 0.0625 within the error bar. 
For the SPACS we obtain a qq a pp = 0.101 ± 0.006 > 
0.0625. The coherent state has the minimal uncertainty 
indeed and, therefore, is pure. This result is in agreement 
with the detection imperfection discussed in Sec. lIIDl bc- 
cause the imperfect detection results in \a) — > \y/fja), i.e. 
the pure coherent state is transformed into another pure 
coherent state which exhibits the same minimal uncer- 
tainty. In fact, this observation confirms the validity of 
using vacuum state for the initial calibration. 

There exists, however, a stronger version of the Heisen- 
berg uncertainty relation which takes into account the 
purity of the state, namely, 



a qq a pp > ft 2 $ 2 (/z)/4, 



(12) 



which is also known as purity-dependent uncertainty rela- 
tion [25|. Here, the purity-dependent function $(/i) = 2— 
V2/x-l if | < fi < 1 and $0) « (4+^16 + 9/?)/9/x± 
4% within the whole range /i <G (0,1]. Employing the 
previously found values of the purity (Sec. UVA"]) . the in- 
equality (JTSJl transforms into 0.101±0.006 > 0.085±0.006 
for the SPACS, which is the first direct experimental 
verification of formula (1121) within the accuracy ~ 3u. 



This result also encourages a feasible verification of two- 
mode uncertainty relations [57| because the correspond- 
ing methods of detecting two-mode states by a single 
homodyne detector are already available |58j |. 



2. State- extended uncertainty relation 

Recently, Trifonov generalized uncertainty relations for 
a pair of different states [26|, [23], where the variances 
of one state were connected with the variances of the 
other by a series of so-called state-extended uncertainty 
relation. One of such relations reads 



V-'qq u pp 



V 2 V (1) 

99 PP 



> 



4 ' 



(13) 

We associate states "1" and "2" with the coherent state 
and the SPACS, respectively. Using the experimental 
data, the relation (Q3) takes the form 0.160 ± 0.006 > 
0.0625, and thus is fulfilled with a great margin. The 
great margin is due to the fact that "1" is a coherent state 
for which a qq = a pp = h/2. This first demonstration 
of state-extended uncertainty relation can encourage its 
further applications to other states saturating it (e.g., 
some squeezed states). 



D. Experimental check of entropic relations 

1. Shannon entropy 

Given a wavefunction ip(q) of some pure state and 
the wavefunction ip(p) of the same state in the mo- 
mentum representation, we are aware that they are not 
independent and are related by the Fourier transform. 
In view of this, the narrower the distribution |V'(9)| 2 
the wider (^(p)! 2 is and vice versa. It means that 
the entropies S q = — f \ip(q)\ 2 \n\ip(q)\ 2 dq and S p = 

— J \ip(p)\ 2 In \ip(p)\ dp cannot take small values simulta- 
neously and turn out to satisfy the following relation [29j : 

S g + S p > ln(wh) + 1, (14) 

which is also valid in case of mixed states, with |-0(g)| 2 
and |f/'(p)| 2 being replaced by the marginal distributions 
w(X, 0) and w(X, 5), respectively. Since quadrature op- 
erators Xg and Xg +v / 2 satisfy the same commutation 
relation as q and p do, one can readily generalize (fi~4|) 
and write f3^-[36Jj 

S(9) + S(9 + 7r/2)>\n(irh) + l, (15) 

where S(6) = Sx e and the right-hand side equals 1.45 
if h = \. The inequality (TT5|) holds true for all LO 
phases 9. Considering 9 as an additional independent 
variable, one can now integrate (fT5|) over 9 £ [0, ir]. Tak- 
ing into account that the theoretical tomogram satisfies 
w(X, 9) = w(—X, 9 + 7r), we obtain 



2 / S{9)— > ln(Trfi) + 1 

Jo 7T 



(16) 



or, equivalently, 



H x ,e = - 



K 



+ OG 



dX 



(If) 



w(X,6)lnw(X,6) 



> i [ln(nh) + 1] , 



(17) 



where Hx,e can be treated as the conventional Shannon 
entropy of the probability distribution function w(X,9) 
of two random variables X e (—00, +00) and 6 £ [0,tt] 
such that i J_ oo dX J d9 w(X, 9) = 1. A similar treat- 
ment of the quantum homodyne tomography as an infor- 
mationally complete positive operator- valued measure on 
[0, 2n] x K is presented in the paper [50|. From the view- 
point of foundations of quantum mechanics, a quantum 
state is defined by a fair probability distribution function 
w(X, 0) of two random variables (a point in the simplex 
of infinite dimension) such that its entropy necessarily 
satisfies the relation (|17p . 

The evaluation of the integral (fTl)]) is performed by 
a trapezoid method, i.e. 2 f* S(0)% « EjlYVj+i - 
°j)[S{0j) + S(9 J+ i)}. The evaluation of S(6j), in its turn, 
is performed by substituting the experimental binned his- 
togram h(X,6j) for w(X,6). 

The finite bin width b is known to affect the right-hand 
side of the relation (fT5"]) (see [601 an d references therein) . 
If we take the bin width b = 0.075 and the cutoff value 
of X equal to 3, then the right-hand side of (TT5|) is to be 
diminished by 0.03 and equals 1.42 for the choice h = i. 
In fact, the allowance is always negative and vanishes 
for larger cutoffs because the states of our interest are 
localized quite close to center of the phase space. Using 
quadratures Xe=o and Xg =n /2, we calculate the left-hand 
side of (fl4|) and the result is 1.43 ± 0.01 for the coherent 
state and 1.65 ±0.03 for the SPACS, where the errors are 
evaluated by comparing the experimental values S{9) and 
S(6 + it). The coherent state saturates the boundary as 
it is predicted by the theory |29| . 

As to integral relation (fT6| , the experimental data yield 
the following quantities of the left hand side of (fT6"|) : 
1.42 ± 0.01 for the coherent state and 1.70 ± 0.03 for the 
SPACS, where the error bars comprise both the error of 
calculation and the errors of the experimental data. 



2. Renyi entropy 

The Renyi entropy of the probability distribution p(X) 
is defined through Rp[p{X)] = (1 - (3)- 1 In J pf l (X)dX 
and represents nothing else but a one-parametric family 
of entropy measures [61]. The Renyi entropy reduces to 
the Shannon entropy in the limit j3 — > 1 . For (3 > 1 there 
exists a conjugate parameter 7 such that fi~ l + j^ 1 = 2. 
We can put j3 = (1 — r) _1 and 7 = (1 + r) _1 , where 
r £ (0, 1). An analog of the relation (Tl"4"|) in terms of the 
Renyi entropy is Rp[w(X,0)] + R^[w(X : f )] > Infrft) - 
|[(1 — /3) _1 In/? + (1 — 7) _1 In 7]. In terms of a single 
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FIG. 6: (Color online) Experimental check of uncertainty re- 
lations for the Renyi entropy (|18|) . Filled circles and squares 
correspond to experimental coherent and SPAC states, re- 
spectively (for points without error bars the corresponding 
errors are greater than 1). Solid line is a theoretical bound. 
When r — > the inequality ()18[) transforms into the inequality 
for the Shannon entropy (|14[) . 



parameter r this relation takes the form [36 



K(r) 



1 + r 



1-r 



hJ [w{X,0)}( 1+rrl dX 

itx^^-^dxX 



111 



> \n(Trh) + i [(1 + r) ln(l + r) - (1 - r) ln(l - r)], (18) 

which remains true by replacing w(X, 0) and w(X, ?) 
by w(X,9) and w(A, -| + 6*), respectively. For the sake 
of simplicity, we concentrate on the experimental check 
of inequality (fl"gj) for r G (—1,1). The results are pre- 
sented in Fig. [SI where the points correspond to the left- 
hand side of (fT8|) calculated via the experimental his- 
tograms. As above, the coherent state saturates the 
boundary (within the experimental errors) and the ex- 
perimentally determined values are symmetrical with re- 
spect to change r — > —r because the position and mo- 
mentum are identically distributed. This does not take 
place for the SPACS and such an asymmetry is readily 
seen. 



V. SUMMARY 

In this paper, we have considered a relatively simple 
but extremely powerful experimental apparatus to mea- 
sure quantum states of light - homodyne detector. Our 
main idea was to use measurable quantities (histograms) 
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to reveal as much information about light states as pos- 
sible. First, the measured histograms enabled us to es- 
timate the optical tomogram, i.e. the quantum state it- 
self. We developed a method for choosing an optimal 
bin width, which ensures statistical confidence and pre- 
vents from undcrsampling at the same time. Second, 
we managed to accomplish a quantitative analysis of the 
accuracy of estimated tomograms by using the peculiar 
property of fair tomograms w(X, 9) — w(—X, 9 + ir). Dis- 
tinction of our approach is that the evaluated errors com- 
prise both statistical and systematical errors. Moreover, 
the detailed analysis can also reveal probable sources of 
systematical errors such as imprecision of the LO phase 
control, which is hardly possible to detect by other meth- 
ods. Even if the systematic error cannot be got rid of, 
one can use an original collection of experimental data to 
postselect those data which exhibit the least systematic 
error. Third, we used the measurable quantities (his- 
tograms) to calculate the characteristics of the state di- 
rectly. For instance, the purity and its error are naturally 
calculated on the basis of measured experimental data 
without any time-consuming state-reconstruction proce- 
dure with controversial error estimation. Our data result 
in the relative error about several percent (l-j-5%) for al- 
most all state characteristics (when their theoretical val- 
ues do not vanish). Last but not least, the operational 
use of the data allowed us to check for the first time the 
fundamental properties of quantum objects such as the 



(purity-dependent) uncertainty relations for position and 
momentum as well as their entropic analogs. 

To conclude, we believe that the developed methods 
will contribute to achieving a higher precision of optical 
homodyne detection and encourage the operational use 
of experimental data, which can turn out to be crucial 
for the analysis of multimode quantum states. 
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